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ABSTRACT 

We study the infrared (IR) properties of galaxies in the cluster MS 1054-03 
at z — 0.83 by combining MIPS 24 /zm data with spectra of more than 400 
galaxies and a very deep K-band selected catalog. 19 IR cluster members are 
selected spectroscopically, and an additional 15 are selected by their photometric 
redshifts. We derive the IR luminosity function of the cluster and find strong 
evolution compared to the similar-mass Coma cluster. The best fitting Schechter 
function gives L* IR = 11.49i jjgA3 with a fixed faint end slope, about one order of 
magnitude larger than that in Coma. The rate of evolution of the IR luminosity 
from Coma to MS 1054-03 is consistent with that found in field galaxies, and 
it suggests that some internal mechanism, e.g., the consumption of the gas fuel, 
is responsible for the general decline of the cosmic star formation rate (SFR) 
in different environments. The mass-normalized integrated SFR within 0.5i?2oo 
in MS 1054-03 also shows evolution compared with other rich clusters at lower 
redshifts, but the trend is less conclusive if the mass selection effect is considered. 
A nonnegligible fraction (13 ±3%) of cluster members, are forming stars actively 
and the overdensity of IR galaxies is about 20 compared to the field. It is unlikely 
that clusters only passively accrete star forming galaxies from the surrounding 
fields and have their star formation quenched quickly afterward; instead, many 
cluster galaxies still have large amounts of gas, and their star formation may be 
enhanced by the interaction with the cluster. 

1 Steward Observatory, University of Arizona, 933 N. Cherry Avenue, Tucson, AZ 85721 
2 NSF Astronomy & Astrophysics Fellow 
3 NOVA Fellow 

4 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

5 Leiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, The Netherlands 

6 National Optical Astronomy Observatory, 950 North Cherry Avenue, Tucson, AZ 85719; Leo Goldberg 
Fellow 



-2- 



Subject headings: galaxies: clusters: individual (MS 1054-03) — galaxies: lumi- 
nosity function — infrared: galaxies 



Introduction 



It has long been known that galaxies in rich clusters are quite different from those in the 



same redshift has been found from local epochs up to z ~ 1 ( 


Kennicutt 


1983; 


Baloe 


l et al. 


1997; 


Hashimoto et al. 


1998; 


Poggianti et al. 


1999; 


Postman et al. 


2001; 


Gomez et al. 


2003 





Several physical mechanisms are proposed to explain this suppression effect, e.g., galaxy- 
gala xy interaction, ram pres sure stripping, and strangulation of the gas reservoir of galaxies 



sec 



Boselli fc Gavazzil 120061 . and references therein). Although all of those processes may 

R in clusters, which one leads to the major effect is still 



play a role in suppressing the S 



controversial (Balogh et al.l 



20021 ; iBaloeh et al. 



2004 



1997 



Has himoto et al.l 1 19981 ; iPoggianti et al.l Il999l ; iLewis et al. 



Kauffmann et al.ll2004l ). More detailed study of the star forming 



activities in rich clusters is necessary to clarify the environmental effects on the SFR. 

Aside from the general suppression compared with the field, the star formation activity of 
galaxies in clusters also evolves with epoch. This behavior was first discovered as the increas- 
ing number of blue galaxies in high z c l uster s compared to local ones , the so-called Butcher- 



Oemler effect (Butcher fc Oemler 



19831 : lEllineson et al. 



2001 



1978 



1984). Spectroscopic studies (iDressler fc Gunnlll982 
Tran et al.ll2005l ) reveal those blue galaxies as star-forming galax- 
ies and indicate increased star forming activity in distant clusters compared to their relatively 
'quiet' local kin. In parallel to this trend in clu sters, the average cosmic star formation rate 
has a l so experienced a rapid decline since z ~ 1 (ILillv et al.lll996t iMadau et al.lll998l ; [Hopkins 



2004: Le Floc'h et al. 2005; Perez- Gonzalez et al 



20051 ) . This raises interesting questions of 



what drives the cosmic decline and how it relates to the decline of the SFR in clusters. 
If clusters only passively accrete star forming galaxies from the field, then the increased 
number of blue galaxies in high z clusters may merely result from the increased number 
of blue galaxies in the field and clusters would show a similar evolution in SFR as in the 
field. Also, if the cosmic decline of the SFR is mainly caused by an overall progressive gas 
consumption in galaxies, the same evolutionary trend in different environments will be a 
natural consequence. On the other hand, if the SFR suppression is largely caused by an 
environment-dependent mechanism, the rate of decline of the cosmic SFR is expected to be 
correlated with the hierarchical growth of structure. 



These questions may lack straight forward answers in the sense that all of the possible 
mechanisms may play a role in shaping the SFR of galaxies we observe, but with importance 



- 3- 



varying with epoch and environment. Howev er, an effort to measure the SFR in clusters at 
different redshifts (e.g.. iPoggianti et al.ll2006l ) and a comparison of clusters with low density 



regions can provide us more insights into these issues. 

Most of the previous studies of star forming activity in rich clusters are based on the 
galaxy colors, UV luminosity and emission-line measurement. The colors of galaxies depend 
on the star formation history, as well as reddening, initial mass function (IMF), age, and 
metallicity. They can only provide an approximate estimate of the star formation history of 
the galaxies. The UV continuum directly traces the emission of the young stellar population 
(< 10 s yr), but it is heavily attenuated by dust and presents only the unobscured star forma- 



tion, which is usually a small fraction of the total star formation in galaxies (e.g.. lBuat et al. 
20051 ). Several emission lines, e.g., Ha, [O II]A3727, provide sensitive, instantaneous mea- 
surement of star formation (SF). The spectroscopic surveys of Ha lines are mostly limited to 
galaxies with z < 0.4, and it is hard t o extend them to a large number of galaxies. The Ha 
narrow band imaging t echnique (e.g.. lKennicuttlll983l ; iGavazzi et al.lll998l ; iBalogh fc Morris 
2000l ; iFinn et al.l I2004T ) is more efficient but it requires specific filters for different redshifts 
and the continuum measurement could bring extra uncertainties. Due to the [O 11}/ Ha 
variations in galaxies and the higher extinction, the [O II]A3727 lines are a less precise SFR 
tracer than Ha, but they can be applied to higher redshift (z > 0.4) galaxies. Both Ha 
and [O II]A3727 emission line measurements are affected by dust extinction, which can be 
very high for strongly star-forming galaxies (e.g.. lChoi et al.ll2006l ). As a result, robust mea- 
surements of the SFR unaffected by extinction in clusters are still lacking, especially at high 
redshift (z > 0.5). 

On the other hand, the infrared (IR) bolometric luminosity from the interstellar dust 
heated by the young stars in ga laxies can provide a sensitive SFR indicator minimally affected 
by extinction (jKennicuttlll998l ) . Moreover, bec ause mid-IR broadband emission shows a good 
correlation with the total IR luminosity (e.g., iTakeuchi et al.ll2005l ; iDale et al.ll2005l ). it has 
become common practice to use single band mid-IR emission as a star formation indicator. 
The mid-IR correlates well over a large lu minosity range with the extin ction corrected optical 



and near infrared SFR indicators (e.g., lAlonso-Herrero et al.ll2006l ). Although there are 



many uncertainties involved in this method, e.g., the escape of the UV photons in optically 
thin regions, the heating of dust due to the ionizing photons from older stellar populations 
and the uncertainties in co nverting single band emission into the total IR luminosity (e.g., 
Perez-Gonzalez et al.ll2006l ). it is still a robust method to measure obscured SF in luminous 
galaxies. 

With data from the European Space Agency's Infrared Space Observatory (ISO), the 
IR-bright galaxies in quite a few clusters have been studied, with many important results 
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(see review in iMetcalfe et al.ll2005l ). In their study of A1689, iDuc et al.l (120021 ) found at 
least 30% of the 15 /im sources show no evidence of current star formation in the optical 
spectra and that at least 90% of the SF in A1689 is obscured by dust. Their study shows 
the importance of using mid-IR emission to detect obscured SF in clusters. However, with 
only a small sample of low to mid-redshift clusters (z < 0.5) studied by ISO, a complete 
understanding of the SF in clusters has not yet been achieved, especially for high redshift 
clusters. 



The Multiband Imaging Photometer (MIPS, iRieke et al.l 120041 ) on the Spitzer Space 
Telescope, with its high resolution and sensitivity, provides efficient measurements of IR 
emission for large samples of galaxies. Using the 24 /im observations by MIPS, we can 
study the SFR of clusters up to z ~ 1. These data provide a good chance to expand and 
explore the study of SF in clusters, leading to systematic and comprehensive understanding 
of galaxy evolut ion in them. We have already studied the IR properties t he Coma cluster 



(iBai et al.l l2006h at z = 0.02 and RXJ0152.7-1357 (RXJ0152 hereafter) flMarcillac et al. 
20071 ) at z = 0.83. Here, we present an IR study of another well known cluster MS 1054-03 
(z = 0.83). At this redshift, MIPS 24 /im data t race the rest frame 13 um flux. The 12-15 ^im 



emission has been demonstr ated particularly by lRoussel et al.l (120011 ) ; lTakeuchi et al.l (120051 ); 
Alonso-Herrero et al.l (120061 ) to be an accurate SF tracer. In §2, we describe the MIPS data 
for the cluster. §3 discusses the analysis of the IR and optical data. We present the IR 
luminosity and additional IR properties of the cluster in §4. In §5, we discuss the results 
and we summarize them in §6. Throughout this paper, we assume a ACDM cosmology with 
parameter set (h,Q ,A ) = (0.7,0.3,0.7). 



DATA 



2.1. Observations 



The MS 1054-03 field was observed at 24 /im by MIPS on June 2005 in photometry 
mode. The total MIPS field has a size of 5' x 10'. The integration time is about 3600 second 
pixel -1 in the central 5' x 5' region, and is about 1200 second pixel -1 in the re st of the field. 



The data were processed with the MIPS instrument team Data Analysis Tool (IGordon et al. 



20051 ; lEngelbracht et al.l 120071 ) and scan- mirror-dependent flats were used. The final mosaic 
has a pixel scale of ~ 1.25" pixel -1 and a point-spread function (PSF) with FWHM ~ 6". 

The spectroscopic data were obtained with the Low Resolution Imaging Spectrograph on 
the Keck Telescope based on the wide-field images taken by the Hubble Spac e Telescope (HST) 
WFPC2 in F606W and F814W Ivan Dokkum et al.llioool : iTran et al.l 120071 ) . They cover the 
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5' x 5' central region of the cluster and yield more than 300 redshifts, adding to a total of 
more than 400 galaxies with previously known redshifts. Among them, 144 sources with 
reliable redshifts are identified as cluster members. The spectroscopic survey is about 50% 
complete down to Jgu = 22 mag and it is mostly coincident with the central deep region of 
the 24 /im observations. 

Along with the HST photometric data, UBV and near-IR J S HK S data in a similar region 
of this field were obtained with FORS1 and ISAAC on the Very Large Telescope (VLT) as 
the part of the Faint InfraRed Extragalactic Survey (FIRES). From these data a K-band 
selected catalog with 1859 sou rces was extracted. The photometric catalog is presented in 



Forster Schreiber et al.l (120061 ) and it is 90% complete to K s> ab ~ 24.1 mag . Photometric 



redshifts (z„ h) were derived from this catalog using the method presented by iRudnick et al. 



(120011 . 120031 ). The photometric redshifts are less accurate than the spectroscopic redshifts 
(z sp ), with Sz/(1 + z sp ) = 0.074 for z < 1, but the FIRES catalog is much deeper than the 
spectroscopic data. Therefore, we use photometric redshifts to supplement the spectroscopic 
catalog. 



2.2. Source Detection and Completeness 



Since the cluster members are not resolved at 24 /im, we used DAOPHOT II (Stetson 



19871 ). a package for PSF fitting photome try, to detect so u rces a nd measure their fluxes. We 
follow the same strategy as described in iPapovich et al.l (120041 ) . Because of the significant 
zodiacal IR emission at low ecliptic latitudes, the 24 fim background level in the MS 1054-03 
field is fairly high, averaged at about 40 MJy sr" 1 . Even with more than 3000 seconds of 
integration time, the detection limit is not as deep as in some low background regions. By 
adding artificial point sources into the image, we found a 80% completeness limit at about 
80 /xJy In the left panel of Fig. [H we plot the completeness of the 24 jum detections vs. the 
flux density. The completeness drops from 80% at about 80 /xJy to only 50% at 50 /iJy. We 
detected about 240 sources with / 2 4 > 50 /xJy in the central region that is covered both by 
IR and optical data, about 180 of them with / 2 4 > 80 /iJy. 



ANALYSIS 



3.1. Spectroscopically Confirmed Cluster Members 



We select the galaxies with 0.81 < z sp < 0.85 as cluste r members, which corresponds to 
a 3-cr line-of-sight velocity dispersion of 1156 ± 82 km s _1 (ITran et al.l 120071 ). We select 144 
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cluster members out of around 400 galaxies with spectroscopic data. Since the spectroscopic 
data are only 80% complete down to Igu = 21, we need to correct for the incompleteness 
to avoid bias. In the right panel of Fig. [U we plot the ratio of galaxies with successful 
spectroscopic redshift measurements among all galaxies in the imaging data as a function of 
isu- We use this curve to correct for the spectroscopic incompleteness when calculating the 
IR luminosity function (LF). 



3.2. Photometrically Identified Candidate Cluster Members 

Even though we can roughly correct for the spectroscopic incompleteness of our IR 
sample using the completeness curve, about one third of the IR luminous members have a 
magnitude of Isia > 22, where the incompleteness is larger than 50% and a simple correction 
can be erroneous. In addition, the igu band selection is biased toward the blue galaxies and 
may miss some dusty star forming galaxies with extreme red colors. This would make the 
incompleteness correction based on the J 8 i 4 magnitudes inadequate. Therefore, we also use 
photometric redshifts to help select the cluster member sample. 

Because most of the cluster members have Isu — K s> ab ~ 1-0, the 90% completeness 
limit (K St AB ~ 24.1) of the photometric survey corresponds to J 8 i4 ~ 25.1, indicated as a 
dashed vertical line in the right panel of Fig. [TJ This limit is about three magnitudes deeper 
than the spectroscopic survey. However, the uncertainties in the photometric redshifts are 
large (8z ~ 0.14 at the cluster redshift z c i = 0.83) compared to those of the spectroscopic 
redshifts and their distribution is non-Gaussian. Because of this, a simple cut in z p % is 
not effective to select cluster members and would cause large contamination. Therefore, we 



use t he probability curve of the z ph deduced from Monte-Carlo simulations (IRudnick et al. 



20031 ) to select the cluster members. If the integrated probability of cluster membership for 
a galaxy with K s> ab < 24.1 over the range of z c i — 0.14 < z < z c i + 0.14 is larger than 60% 
(normalized by the total probability), we designate this galaxy as a cluster member. The 60% 
threshold is selected to best balance between maximizing correct selections and meanwhile 
minimizing incorrect selections when applied to galaxies with spectroscopic measurements. 
When the threshold is set to 70%, the incorrect selection drops from 17% to 15%, but the 
correct selection also drops from 83% to 73%. Altogether, we select 454 candidate cluster 
members from 1858 sources in this region, three times more than the members selected by 
spectroscopic redshifts alone. 

The performance of the photo-z selection is expected to decrease for faint K band 
sources. As the uncertainties in the fluxes increase, the internal uncertainties in the photo-z 
increase too. Therefore, the probability curve of z p h broadens and its integrated value in the 
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same redshift range will be lower. As a result, more faint galaxies will be rejected given the 
same threshold. 



3.3. Crossmatch between Optical/NIR Sources and 24 /mi Sources 

In a crowded field such as MS 1054-03, we need to be careful in crossmatching between 
the Optical/NIR sources and the 24 /im sources. The optical and N IR images have an 
absolute astrometric accuracy of < 0"5 iForster Schreiber et al.l booeh. The astro metry 



of the 24 /zm image is calibrated using the USNO-B1.0 catalog ( Monet et al.l I2OO3I ) . and 
has an accuracy of rms < 0'.'6. However, due to the rather large FWHM (~ 6") of the 
24 /iin image, this accuracy in position can only be achieved for bright IR point sources, 
namely with a 24 /im flux > 100 /iJy. For faint sources, as shown by our simulations, the 
average uncertainties of the positions are about 1". We use a radius of 2" (~ 15 kpc at 
z = 0.83) to correlate the optical/NIR cluster members with their IR counterparts. This 
matching radius accounts for the possible displacement between the optical/NIR and 24 
/im brightness centroids, the astrometric uncertainties and local astrometric offsets. If more 
than one counterpart is found in this radius, the nearest one is picked. We estimate the 
chance of random matches by randomly re-distributing the IR sources and matching them 
with the same criteria. We found only 4.8 ± 0.9% random matches. Using this criterion, we 
obtain the preliminary matching lists for both spectroscopic and photometric samples. We 
then carefully check each individual source by eye to exclude any apparent mis-identification, 
e.g., contamination from nearby bright IR sources. Finally, we identify 19 sources selected 
by spectroscopic redshift and an additional 15 sources selected by photometric redshift with 
IR emission > 50 /iJy. We refer to those 19 IR galaxies with spectroscopic redshifts as our 
spectroscopic sample, and those 19 sources plus 15 sources selected by photometric redshifts 
as our combined sample. For the 15 photometrically selected galaxies, we will use the cluster 
redshift as their redshifts in the following study. In addition, for most of these IR galaxies, 
we use the IRAC data kindly provided by the FIRES group to confirm the crossmatching. 



3.4. Incompleteness Correction 

Many of the sources in our final samples are fainter than the 80% completeness limit of 
the spectroscopic and 24 /im surveys. It is therefore necessary to correct for the incomplete- 
ness of both surveys to have an unbiased number count. To do this, we use the inverse of the 
completeness curves in the 24 /im and / 814 bands given in Fig. [T] as the weighting functions 
to calculate LFs. For the 19 galaxies in the spectroscopic sample, we correct both for the 
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spectroscopic incompleteness and 24 /im incompleteness according to each galaxy's Igu mag- 
nitude and 24 /im flux density. All the galaxies in the combined sample are brighter than the 
90% completeness limit of the photometric survey, so we only correct for the incompleteness 
in the 24 /im detections. 

The incompleteness correction can be very large (see Fig. [1]), especially for faint galaxies 
in the spectroscopic sample. It boosts the number density up to 3 times at the faint end of 
the luminosity function. We will discuss the effect of the incompleteness correction in §4.1. 



3.5. Deduction of the Total IR Luminosity 



To maintain continuity with iBai et al.l (120061). we use the i r met hod to determine total 



IR luminosities. We shift the SEDs given by iDevriendt et al.l (119991 ) to the cluster redshift 



to deduce the rest-frame total IR luminosities (Lj R , A = 8 — 1000 /im) of the galaxies. Those 
SEDs are based on a sample of nearby galaxies and include three types: normal spirals, 
luminous IR galaxies (LIRGs) and ultraluminous IR galaxies (ULIRGs). The deduction 
of the total IR luminosity depends primarily on the ratio between the rest frame Lir and 
the 13 /im luminosity (£24/(1+2) ~ £13) for each galaxy. The template SEDs indicate that 
this ratio is almost constant within each type, but increases by three times from normal 
spirals to ULIRGs (see Fig. [2]). Since the K s band (similar to rest frame J band) flux and 
24 /im (similar to rest frame 13 /xm) flux are good indicators of the old and star-forming 
components of galaxies respectively, the color between these two bands (//^//W) can be 
used to distinguish different types of galaxies. In Fig. [21 we show the correlation between 
Lir/L 2 4 : /(i +z ) and JkJ fi<± for each type of SED. The open stars are ULIRGs, the square is 
a LIRG, and the open triangles are normal spirals. We interpolate the fx s /f2A colors of the 
cluster members linearly into the correlation given by template SEDs and get a Lir/Lh/q+z) 
ratio for each galaxy. The filled circles are the spectroscopic sample and the filled upside 
down triangles are the combined sample. According to their fK s / f 24, colors, about one third 
of the galaxies in our spectroscopic sample and about half in the combined sample have 
ULIRG or LIRG SEDs. Using the Lj#/L 2 4/(i+z) ratio given by the interpolation, we deduce 
the total IR luminosity from the 24 /im flux of each galaxy in our samples (see Table 1). 

The method we used above basically assumes that there is no intrinsic variation in SEDs 
for galaxies with same fic a /fM colors and that the templates represent a complete sample 
of IR galaxies up to z ~ 0.8. However, both of these assumptions are questionable, given 
the large variation of IR SEDs among star forming galaxies and the possible evolution of 
galaxy properties from z = 0.8 to z — 0. To estimate the uncertainties of Ljr caused by the 



limitations of the SED templates, we used a different set of SEDs from lDale &: Heloul (120021 ) 
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and a strategy descri bed in iMarcillac et al.l (l2006bl ) to deduce the total IR luminosities. 



Marcillac et al.l (120071 ) use this method to deduce the total IR luminosities in another cluster, 
RXJ0152, at a similar redshift. We plot the deduced L IR / '£ 2 4/(i+ z ) ratio of galaxies vs. 
fx a l f-2A colors as small dots in Fig. [2] for comparison, though the analysis itself does not 
depend on fxs/fii color. This method gives a slightly smaller typical Lir compared to the 
first method, by a factor of 0.9 ±0.3 on average. The difference is more pronounced for those 
galaxies with a smaller fx a / f24 color (< —1), where the difference is up to a factor of 2, 
and may be caused by the mis-classification of SED types with only one color. To exclude 
this possibility, we compared the multi-wavelength photometr y of sample galax i es (op tical 
+ NIR + IRAC + MIPS 24 /im ) with the model SEDs from bevriendt etaD Jl999h and 
confirmed they do have LIRG/ULIRG type SEDs. The difference, caused by the wide SED 
variations from galaxy to ga laxy, is typical of methods to estimate to tal IR luminosities from 



24 /im measurements (e.g., iPapovich fc Belli 120021 ; iDale et al.l 120051 ) . It does not affect the 



results of this paper significantly. If we do not consider the uncertainties caused by the SED 
fitting, the error of Lm is dominated by the flux uncertainties at 24 /im, which are typically 
50% for the galaxies studied in this paper. 



3.6. Contamination from AGNs 



When we deduce total IR luminosities for the galaxies, we assume their IR emission 
is entirely from emission by dust heated by star forming activity, neglecting the possible 
contribution from active galactic nuclei (AGNs). Although optical studie s suggest that 
AGN s reside in only about one percent of galaxies in clusters up to z ~ 0.5 (IDressler et al. 



19991 ). recent X-ray surveys have found an excess of point sourc es in cluster fields, many of 
which are confirmed as cluster AGNs (e.g.. iMartini et al.l 120021 ) . These discoveries suggest 
that AGN contamination may be an issue, especially for our small number samples. 

In the MS 1054-03 field, s urveys in the rad io and the X-ray bands have been analyzed 
to identify the possible AGNs. iBest et al.l (120021 ) conducted an extremely deep 5-GHz radio 
observation and found 34 radio sources, 8 of which are confirmed as cluster members by 
their spectroscopic data. On the basis of the [O II] emission line flux and radio flux density 
ratio, they further conclude that 6 of these 8 radio sources are AGNs, on e source (No 5) is a 
star-forming galaxy and one (No 14) is ambiguous. iJohnson et al.l (120031 ) analyzed the 91 ks 
Chandra observations of the cluster and detected 47 X-ray sources. Among them, two sources 
are confirmed as AGNs (source 7 and source 19); source 19 is also detected in the radio. 
Altogether, there are 8 possible AGN members in the cluster. To avoid losing possible AGN 
candidates in our IR galaxy sample due to the incompleteness of the spectroscopic survey, 
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we crossmatched the IR galaxies with all the 34 radio sources and the 47 X-ray sources. 
The two confirmed AGNs (X-ray source 7 and source 19) and the No 5 radio source are 
detected in the IR. We exclude the two AGNs from both of our sampl es. For the No 5 radio 
source, we use a radio spectral index of -0.8 and the formula given by lHopkins et al.1 (120031 ) 
to convert the radio flux to the SFR. We deduce a SFR of about 88 M Q yr" 1 from its radio 
flux, which is consistent with the SFR e stimated from the total IR luminosity, ~ 61 M yr _1 , 
using the conversion formula given by iKennicuttl (119981 ) . This agreement further confirms 
radio source No 5 as a star-forming galaxy. Due to the limitations of the spectroscopic data 
and the sensitivity of the X-ray survey, we can not totally exclude all AGN contamination 
from our IR galaxy samples, but the fact that only one star forming IR galaxy in our sample 
is detected either in the radio or in the X-ray band indicates the contamination is negligible. 



3.7. Comparison of the IR and [O II] emission line deduced SFRs 



Many previous studies of the SFR in clusters at z > 0.4 rely on the [O II] emission line 
as an indicator, which is very sensitive to extinction and metallicity. We compare the SFR 
deduced from the [O II] emission line luminosity and that from the total IR luminosity in 
Fig. [3] for the IR galaxies in the spectroscopic sample. Among 15 IR galaxies with [O II] 
data, 12 galaxies have emission lines. The [O II] emission line luminosity is estimated by 
multiplying the equivalent width of the emission line by the continuum flux. The continuum 
flux at the rest frame of the [O II] line is approximated by the continuum flux in Vqo6- The 
[O II] emission line luminosit y is converted to a SFR using the formula SFR[ n] = (6.58 ± 



1.65) x 10 42 L [on] (ergs s *) (IKewley et al.1 120041 ). where L[on] is the luminosity corrected 



for extinction. Without any extinction correction, the SFR[qh] is smaller than the SFRir 
by more than one dex on average (the open circles), but with large scatter. Since we do 
not have enough optical data to deduce extinction, w e used the I R-lum inosity-dependent 



extinction A l y = O.751og(LiR/L ) — 6.35 mag given by IChoi et al.l (120061 ) to correct for the 



dust attenuation. This extinction formula is deduced from the ratio of the SFRjr and the 
SFR measured from emission lines, assuming SFRjr approximate the true SFR. The galaxies 
in our sample all have a Ay 1 greater than 1.5. We convert A l y to the exti nction of [O II] lin e 
using the reddening curve of ICalzetti et al.l (120001 ). the same one used by lChoi et al.1 (120061 ). 
The extinction-corrected SFR[ n] agrees with the SFRir reasonably well (the filled circles), 
with a scatter of about 0.5 dex. 

For the three IR galaxies without emission lines, we plot the SFRir as the upper limits 
of their SFRroni (the open circles with downward arrows). There are also about a dozen 
[O II] emission line galaxies (EWjoiij > 5A) not detected at 24 /im. For these galaxies, the 



- 11 - 



lack of IR emission probably suggests relatively less dust and smaller extinction, so we used a 
fixed Ay = 1.0 to deduce their SFRpii]- Their SFR[on] are all quite small, with a maximum 
value of 7 M yr _1 , well below the 80% detection completeness limit. 

We also compar ed the SFRpRj of galaxies in RXJ0152 with the SFRjon] given by 



Homeier et al.l (120051 ). correcting for extinction with A l y-. The open and filled triangles 
in Fig. [3] are the data without and with extinction correction. Again, the extinction cor- 
rected SFR[qtti shows a better consistency with the SFR IR . This agrees with the results of 



Marcillac et al.l (120071 ). who also found a large amount of dust-embedded SF in RXJ0152. 



The median values of Ay 1 for the IR galaxies with [O II] emission lines in MS 1054-03 and 
in RXJ0152 are both about 2, corresponding to a correction factor of ~ 14 for the SFR[on]. 
This result implies the star-forming galaxies in these clusters are enshrouded heavily by dust 
and the SFRpii] without extinction correction only measures a small portion of the total 
SFR. Even with the widely adopted extinction of 1 mag at Ha (Ay ~ 1-2), the emission line 
SFRs still underestimate the SFR by a factor of 4 for these IR bright galaxies. 



3.8. Comparison of the IR and Ultraviolet continuum deduced SFRs 



UV luminosity is also widely used to estimate the SFR of galaxies. Although it is 
very sensitive to the dust extinction, it gives us access to the "non dusty" star formation 
and is therefore complementary to the IR-deduced SFR. For the cluster members, we have 



;/,2200A 



from the galactic extinction corrected 



derived the rest-frame 2 200A luminosity (L 

U — K photomet r y from iForster Schreiber et al.l (120061 ). using the methodology presented 
in Rudnick et al. (2003). We estimate the SFR from L ° using the formula given by 



Kennicuttl ( 119981 ). The conversion assumes a Salpeter IMF and a constant SFR, with UV 
emission dominated by a stellar population younger than 100 Myr. These assumptions are 
consistent with t hose used to deduce the IR SFR conversion formula we adopt in this paper 
jKennicuttJfl998h . 



In Fig. HI we show the comparison between the UV continuum-deduced SFRs and the IR- 
deduced SFRs. For the IR bright cluster members, the unobscured star formation are only a 
small fraction of the total star formation. The median value of SFRjr/SFRuv is ~ 12 for the 
spectroscopic sample and ~ 16 for the combined sample. Such a large SFR IR /SFRuv ratio is 
mainly due to the fact that our IR data is only sensitive to galaxies with SFR > 10 M & yr _1 , 
where extinction is known to be large. The combined sample has a larger SFRir/SFRuv 
ratio on average because the spectroscopic survey is J-band magnitude limited and is biased 
against the most dusty star forming galaxies. If we calculate the visual gas medium extinction 
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from Ljr using the formula given by I Choi et al. (all the extinctions we mentioned in 

§3.7 are for the gas medium) and apply the corresponding UV stellar continuum extinction 
to SFRuv, we will have a better agreement between SFR IR and SFRuv, as shown in Fig. HI 
However, even after this extinction correction, there are still many galaxies, especially the 
ones selected by photometric redshifts, showing a much smaller SFRuv compared to SFRtr. 

If we assume SFRtr is the total SFR and directly estimate the NUV extinction by 
2.51og(SFR IR /SFRuv), we will have a median A NUV of 2.7±0.5 mag and 3.0±1.0 mag for the 
spectroscopic sample and the combined sample respectively. The high extinction we found 
here supports our assumption that SFRir provides a reasonable estimate of the total SFR. 
The NUV extinction of the stellar continuum c an be translated into the visual extinction 
of the gas medium using the reddening curve of ICalzetti et al.l (120001). Ay = 2 .9 ± 0.5 mag 
and A v = 3.3 ±1.1 mag respectively for our two samples. iBuat et al.l (120071 ) studied the 
extinction of a sample of LIRGs detected in the Chandara Deep Field South at z = 0.7 
using the ratio of the total IR and FUV luminosity. They found an average FUV extinction 
of 3.33 ± 0.08 mag for their sample, corresponding to a gas medium visual extinction of 
Ay = 2.97 ± 0.07 mag0 , which is in very good agreement with our results. 

In addition to the IR cluster members, we also calculate the SFRuv for all the cluster 
members (spectroscopic + photometric) without detectable IR emission. Their SFRuv are 
all at least two times smaller the detection limit of the SFRtr, which confirms there is no 
galaxy with a high level of star formation that is missed by IR selection due to the lack of 
dust. 



RESULTS 



4.1. IR Luminosity Function 



After obtaining the total IR luminosity of each galaxy, we calculate the LF for each 
sample. For the spectroscopic sample, we correct the number counts for the incompleteness 
in both the igu and 24 /im bands. We only correct for incompleteness in the 24 /zm detections 
for the combined sample. The overlapping area between the spectroscopic survey and the 



lr To be consistent with Choi et al. ( 200(jh . we use the reddening curve of Calzetti et al. (2000) for all 



the extinction conve rsion in this paper , which indicates Ay — 4.05E(B — V) g . It is different from Ay — 



3.1E(B—V) g used bv lBuat et al.l (|2007l ) . If we adopt their conversion method, the absolute value of the visual 
extinctions corresponding to the NUV and FUV extinctions will change, but the results of the comparison 
will remain the same. 
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24 fim observations is about 4.8 Mpc 2 , and the overlapping area between the photometric 
survey and the 24 /im observations is 5.5 Mpc 2 . 

The IR LFs are shown in Fig. The open circles are the LF deduced from the spec- 
troscopic sample with out any incom pleteness correction and the error bars are estimated 



by Poisson statistics (jGehrelsl Il986l ); the filled circles are the results corrected for the in- 
completenesses in both the spectroscopic and IR surveys. The correction is quite significant 
except for the brightest data point. The error for the incompleteness corrected LF is ob- 
tained by multiplying the original error by the incompleteness correction made at each data 
point. Since we do not consider the error caused by the incompleteness estimate itself, the 
error bars should be considered to be lower limits to the actual errors. Similarly, the incom- 
pleteness uncorrected and corrected data points of the combined sample are shown as open 
and filled squares respectively. Even though the uncorrected LFs of the spectroscopic and 
combined samples exhibit a large difference, their incompleteness corrected ones agree with 
each other quite well. This good agreement demonstrates that neither the simplified spec- 
troscopic incompleteness correction nor the uncertainty of the photometric redshifts affects 
our resulting IR LFs significantly. It also shows that there are few galaxies with extremely 
red optical-IR colors missed by the selection limit in the 7 8 i4 band. 

Despite this general agreement, the difference in the brightest data point may cause 
quite a large discrepancy when we try to fit the LF. It also raises questions about the incom- 
pleteness correction because we expect it to be least significant for the brightest galaxies. 
Two spectroscopically and two photometrically selected IR galaxies contribute to this data 
point. The two galaxies selected by spectroscopic redshifts are both very bright (/gu ~ 20) 
late type galaxies, and the two selected by photometric redshifts are both about two mag- 
nitudes fainter in the igu band and slightly brighter in the IR. Even though the probability 
of the photometric redshifts of those two galaxies falling into the one sigma error range of 
the cluster redshift is more than 70%, their best fitting photometric redshifts are both about 
0.95. As an independent check, P. G. Perez- Gonzalez helped us get an other set of photomet- 



ric re dshifts for these two IR galaxies using a different fitting strategy (IPerez- Gonzalez et al. 



20051 ) and with IRAC photometric data as an addition. These photometric redshifts have 
an average accuracy of Az = 0.08. The best fitting redshifts of those two galaxies are 
0.97 ± 0.09 and 1.00 ±0.11, both more than 1 o above the cluster redshift. Their extremely 
large SFRpr/SFRuv ratio, as shown in Fig. HI also suggest them as background sources. 
In addition, among the 20% of photometric sources with spectroscopic data, one galaxy as 
bright as those two at 24 /im is selected as a cluster member by its photometric redshift but 
shown to be a non- member by its spectroscopic redshift. Statistically, it is possible that four 
more foreground or background contaminations may occur in the whole sample. However, 
we still can not rule out the possibility of those two sources as real cluster members given 
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the uncertainties of the photometric redshifts. Spectroscopic data are needed to clarify the 
ambiguity. 

Because the incompleteness-corrected IR LF of the spectroscopic and combined sample 
are generally consistent, while the combined sample has larger uncertainties at the brightest 
end (the spectroscopic sample, on the other hand, should have the smallest uncertainties due 
to the incompleteness correction at this point), we select the incompleteness-corrected IR 
LF of the spectroscopic s ample as the IR LF of the cluster. We fit this LF with a Schechter 
function (]Schechterlll976l ) . Since we only have thr ee data points, w e fix the faint end slope to 
the same value as the IR LF of the Coma cluster (IBai et al.ll2006l ). We adopted a chi-square 
minimization method for the fitting. We also use the non- detection of the brighter galaxies 
beyond the brightest bin as a constraint during the fitting (IBai et al.ll2006l ). The best fitted 
parameters are: 

a = 1.41 (fixed); \og(L} R /L G ) = 11.491°;*°. (1) 

The resulting fit is shown as the solid curve in Fig. [51 Fitting the incompleteness-corrected IR 
LF of the combined sample gives an even larger L* IR , with \og(L* IR /L Q ) = 11.7318;^. We only 
use the Poisson statistical errors for the fitting and do not consider the uncertainties caused 
by the errors in Lj R estimation. The best-fitting parameters have large uncertainties, because 
the Schechter function fitting depends strongly on the brightest bin, which only includes two 
galaxies to constrain L* IR , and even small changes of the Lm of those galaxies can cause large 
changes in the best fitting parameters. We also note that due to the degeneracy between 
the faint end slope and the characteristic IR luminosity (L* IR ), the best-fitting L* IR value we 
obtained here depends on the assumed faint end slope. If we vary the faint end slope from 
its current value by ±0.2, the best-fitting L* IR would vary by ±0.13. However, by fixing the 
faint end slope and fitting to a Schechter function, we can quantify the difference between 
LFs. Because no available infrared data in this redshift range penetrate significantly below 
the LIRG range, virtually all studie s use a fixed low lumi nosity slope, including the field IR 
LF we compare with in this paper (ILe Floc'h et al.ll2005l ). 



We have tested the dependence of these fits on the uncertainties in the deduction of 
L IR with different methods. If we use the L IR of the spectroscopic sample deduced from the 
second method listed in Sect. 3.5, the IR LF does not change significantly, as shown by open 
and filled triangles in Fig.[5J The best fitting function has a smaller \og(L* IR / L & ) = ll^l^;^. 
The difference is still within the one sigma Poisson error, suggesting that small number 
statistics dominate the uncertainty to define a best fitting LF, and the systematic error 
caused by different L IR deduction methods is negligible. Therefore, in the rest of the paper, 
we only use the results from the first method. 
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4.2. Comparison with Coma IR LF 



We compare the I R LF of MS 1054-03 t o that of the Com a cluster, which has similar 
mass as MS 1054-03 (ILokas fc Mamonl 120031 ; iJee et al.l 12005a] ) and whose galaxy infrared 
luminosities are deduced using the same set of SEDs as the first method in this paper. We 
plot the best fitted Schechter function of the Coma cluster IR LF as the dotted curve in 
Fig. [5j The characteristic IR luminosity in MS 1054-03 is ten times larger than that of 
the Coma cluster (\og(L* IR / L Q ) = 10. 49^0:2!) • The surface density of the IR galaxies with 
logLm > 43.5 expected from the MS 1054-03 LF is about 5 times larger than that in Coma. 
We integrate the best fitted Schechter function of the MS 1054-03 IR LF from logL/^ = 44 
to logLjR = 46 and get a SFR density of 190 M & yr" 1 Mpc~ 2 , about 16 times larger than 
the SFR density of the Coma cluster (~ 11.4 M yr _1 Mpc~ 2 ). 

The significant difference between the IR LF s of MS 1054-03 and of Coma agrees with the 



gener a l evolution trend found in the field IR LF (ILe Floc'h et al.ll2005l ; iPerez- Gonzalez et al. 



2005h . ILe Floc'h etaJ (b()05h quantified the evolution of the IR LF in the CDF-S field in both 



density and luminosity as [L* IR oc (l + z) aL , <p* IR cx (l + z) ao ], with the best fitting parameters 



«l = 3.15 ± 1.6, an = 1.02 ± 1.6. Corresponding to this, we estimate the difference of the 
two cluster IR LFs using the same parameters and get «l = 4.0^12, «d = 1-4. For «d, we 



do not give an error estimate due to the large uncertainties of the best fitting (p* IR value. 
The results agree within the errors and even suggest a slightly stronger evolution of these 
two cluster IR LFs compared to the field IR LF. This is also demonstrated in Fig. [51 where 
the dashed curve corresponds to the LF of the Coma cluster evolved to z — 0.83 using the 
field IR LF evolution law. Both the incompleteness-corrected IR LFs of the spectroscopic 
and combined sample of MS 1054-03 fall above this curve. 

The agreement between the evolution of the IR LF in these two clusters and the evolution 
in the field might suggest that the population of star forming galaxies in clusters is dominated 
by recently accreted field galaxies. Therefore, their IR LFs would not be very different from 
that of field galaxies and they would show similar evolution. However, this explanation is 
not favored by our following analysis (see §4.3). More likely, the similarity in the evolution 
trend suggests that the cosmic SFR decline is caused by some general mechanism existing 
both in cluster and field environments, probably the consumption of the gas fuel for SF. 
Nevertheless, due to the intrinsic variation of cluster properties and with only two clusters 
in this comparison, these results are far from conclusive. More clusters need to be studied 
to confirm the evolution trend further. 
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4.3. Comparison with Field IR LF 



The IR LFs we deduced for the cluster are all projected LFs. If we assume the cluster 
has a radius of about 5i?2oo (-R200 is the radius within which the mean cluster density is 
200 times the critical density of the universe at that red shift [f ), we can calcu late the IR 
LF per volume and compare it with the field IR LF from lLe Floc'h et all (120051 ) at similar 
redshift. 5i?2oo is near the turnaround radius of the cluster, and there will be few infalling 
galaxies beyond it. Although the spectroscopic sample is almost free of contamination from 
field galaxies, the redshift selection (0.81 < z < 0.85) still could include a few foreground 
and background field galaxies whose redshifts fall within the cluster velocity dispersion. To 
exclude this field contamination, before we convert the projected cluster IR LF to LF per 
volume, we calculate the projected field IR LF in a cylinder with a length corresponding to 
z = 0.81 to z — 0.85 and subtract it from the projected IR LF of the cluster. In Fig. [6l the 
filled circles are the incompleteness-corrected IR LF per volume of the spectroscopic sample 
after the field subtraction. If we integrate the cluster LF in the range of 10 10 ' 8 < Lj R < 
lO 12 !/^, it shows an oyerden sity of about 21 compared with the IR LF from the CDF-S field 
from iLe Floc'h et al.1 (120051 ). The overdensity calculated here can be affected by the cosmic 
variance from field to field, especially the CDF-S field, i n which a lower g alaxy density up 
to a factor of two is found compared with other fields (jWolf et al.l 120031 ) . However, as a 
first-order correction, ILe Floc'h et al.1 (120051 ) already normalized their IR LFs by the ratio 
between the S-band luminosity densities in the CDF-S and over the t hree fields of COMBO 



17 (Classifying Objects by Medium-Band Observations in 17 filters; IWolf et al.l 12003 . 12004 
Bell et al.ll20o3 ). 



Although the estimate of the actual value of the overdensity has some uncertainties, it 
is clear that there is an excess of IR galaxies in the cl uster compared with the field. Such 
an excess of MIR sources is also found in RXJ0152 (IMarcillac et al.l 120071 ). as well as in 
two intermediate r edshift clusters CI 0024+16 and MS 0451-03 (with a smaller significance) 
(jGeach et al.1 120061 ) . Although the cluster shows a clear overdensity of the IR galaxies com- 
pared with the field, it is still smaller than the overdensity of the cluster in the optical bands. 
A fairer test to examine the star formation level in different environments is to compare the 
fraction of IR galaxies in the cluster and in the field. Among the 144 spectroscopically con- 
firmed cluster members, 19 have 24 jum emission brighter than 50 /xJy and 6 of them have 
L IR > 1O 11 L . There are two AGNs also with 24 /im emission, but we already excluded 
them from the sample. Therefore, the fraction of star-forming galaxies with / 2 4 > 50 /xJy 



2 -R200 is widely used as an approximation to R V i r iai- Strictly speaking, R V i r iai i s closer to i?i nn ■ The ratio 
of R2Qo/Rviriai depends on the mass distribution of a cluster. For a NFW profile (jNavarro et al.lll997l ) with 

RlOO/fs = 6, R2OO /Rvirial ~ 0.75. 
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in the cluster is about 13 ± 3% and the fraction for LIRGs is 4 ± 2%. These fractions are 
barely affected by the incompleteness of the spectroscopic survey. There are 211 cluster 
member candidates selected by the photometric redshifts with K s ab < 22 (approximately 
the detection limit of the spectroscopic survey of I S u = 23) and 15 ± 3% of them are IR 
bright and 5 ± 2% are LIRGs. The results are therefore consistent with the fractions based 
only on the spectroscopically selected cluster members. For the CDF-S field, we select the 
galaxies with 0.81 < z ph < 0.85 and R < 22 .6 mag using the photometric redshifts given 
by the COMBO-17 survey (jWolf et all 12004 ) . The cut in R magnitude approximates the 
detection limit of the spectroscopic survey in the cluster, J 8 i 4 = 23. Altogether, we select 62 
galaxies in an area of 775 arcmin 2 (~ 20 times larger than the cluster field) and 39 of them 
have /24 > 50 /iJy. Two of those field IR galaxies are classified as QSOs. So, the fraction of 
the star forming galaxies in the field sample at the same 24 /zm threshold and of a similar 
redshift range is about 60 ± 12%, much higher than the fraction we found in the cluster. 
Although the photometric redshifts we used to select the field sample are less accurate than 
spectroscopic redshifts, the fraction of the star forming galaxies in the field does not change 
significantly in a large redshift range (0.7 < z p h < 0.95) and therefore redshift uncertainties 
will have little effect on the comparison. The smaller fraction of IR bright galaxies in the 
cluster compar ed with the field is consistent with t he results given by studying the galaxy 
emission lines (ILewis et al.ll2002l ; iGomez et al.ll2003l ) and suggest that galaxies in the cluster 
have a lower level of star formation on average. 

Even though the cluster galaxies have a lower level of SF on average compared with the 
field, a fraction of 13±3% of star forming galaxies is still very substantial considering the short 
IR bright phase, especially for the 6 LIRGs which c onstitute 4±2 % of all the spectroscopically 
selected cluster members. iMarcillac et all (j2006al ) analyzed Balmer absorption lines and the 
4000A break of a sample of LIRGs at z ~ 0.7 and found that the duration of the LIRG phase 
is most likely ~ O.llgjg Gyr. The timescale of star formation in local IR bright galaxies 
is even smaller, ~ 10 7 yr (IGao fc Solomon! 12004 ). If these active galaxies are due solely to 
infall from the field, such a short timescale would mean the cluster would have to accrete 
about 60ljj? LIRGs from the field per Gyr to sustain the observed SF level. This is about 
half of the current spectroscopically selected cluster sample. Even if we consider that abou t 
30% of LIRGs are experiencing their second star bursts in 1 Gyr ( IMarcillac et al.ll2006al ). 
it would still mean more than half of the cluster members are the LIRGs accreted from 
the field in the last Gyr. Such a large accretion rate is very unlikel y . The average smooth 
growth of the cluster masses from the simulation of iRowley et all (12004 ) in the one Gyr 
period (z ~ 1.1 — 0.83) ranges from 10% to 40%. 



However, it is possible that we are seeing a large fraction of star forming galaxies in this 
cluster due to a temporary rise of the accretion rate caused by a major merger/infall event. 
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The quadrupole-like temperature structure and the lack of shock-heated regions between 
the two X-ray peaks of this cluster suggest that the major clumps (the central and western 
clumps) are probably at a postmerger stage and the lack of an X-ray peak in the eastern 
clump may also suggest a recent infalling/passing-by (IJee et al.ll2005al ). These merger and 
infalling events might have introduced a large number of field galaxies into the cluster in a 
short time and boosted the accretion rate temporarily. However, the spatial distribution of 
the IR galaxies does not seem to support this scenario. Although MS 1054-03 has several 
subclumps and one of them has relatively enhanced SF, most of the IR galaxies do not 
concentrate in subclumps; instead, they tend to scatter around the cluster and avoid the two 
major clumps (see §. 4.4). Similarly, the cluster members with 15 /j,m emissions in CI 0024, 
a mid-redshift clus ter with high star formation level, also do not show spatial concentrations 
(jCoia et al.ll2005al ). Such a spatial distribution does not support a major merger /infalling 
event. In addition, RXJ0152 also has a si milar fraction of IR b right galaxies, and they show 
no sign of concentration into a subclump (IMarcillac et al.ll2007l ). This evidence suggests that 
such high fractions of IR galaxies in high redshift clusters might be quite common, and it 
is unlikely that they are all due to major infall events. More likely, these IR galaxies have 
been in the cluster for quite a while. However, they probably have never been close to the 
high density region before and still retain a large amount of gas. The recent SF in these 
galaxies may be triggered either by interaction with the cluster intergalactic medium (IGM), 
with other galaxies, or by tides. In support of this hypothesis, evidence h as been found 
previo u sly for star fo r ming burst s in infalling galaxies into clusters by, e.g., iGavazzi et al. 



fl2003fUCortese et all (120061 ) and Mercurio et all fl2004h . 



An alternative possibil ity is that the L I RGs in the cluster have lasted much longer than 
the time scale estimated by IMarcillac et al.l (j2006al ). If the timescale is an order of magnitude 
longer (~ 1 Gyr), then the accretion rate would be 10 times smaller and would not raise the 
problem of too rapid growth. However, such a long timescale would indicate the accreted 
galaxies could retain their gas for a long period and keep their star formation untouched by 
the cluster environment. This alternative view is again inconsistent with the passive scenario 
that the star formation of field galaxies is quenched quickly after they are accreted into a 
cluster. 



4.4. Spatial Distribution of the IR Galaxies 



The spatial distribution of the IR galaxies in the cluster may help us understand the 
effect of the cluster environment on the galaxy SFR. As indicated by the X-ray and optical 
light distributions, the morphology of MS 1054-03 is quite complex. IJee et al.l (j2005al ) recon- 
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struct a high-resolution mass map of the cluster through ACS weak-lensing analysis. They 
confirm the three dominant mass clumps in the c luster previously reported by a WFPC2- 
based weak-lensing analysis (IHoekstra et al.ll2000l ) and find some detailed substructures for 
the first time. In Fig. [3, we overplot the IR galaxies on this mass contour map. The mass 
map is constructed in units of the dimensionless mass density k, and k > 0.1 corresponds to 
a significance of > 3 sigma. The squares are from the spectroscopic sample, and the triangles 
are those additional members selected by photometric redshifts. The sizes of the symbols are 
proportional to the IR luminosities. For clarity, the LIRGs are also indicated by black dots. 
The three major clumps, eastern, central and weste rn (E, C Sz W), as well as the four minor 
clumps (M1-M4), are labeled on the plot following |Jee et al.l (j2005al ). One distinct feature 
of Fig. [7] is that many IR galaxies are distributed in the outskirt region of the cluster. Two 
thirds of the IR galaxies are located in the low density region with k < 0.1, and this ratio 
could be higher if we take into account projection effects. 

Another interesting feature of the distribution is the lack of IR galaxies in the western 
clump and the southern extension compared to the rest of the major structure. We divide 
the major body of the cluster into two approximately equal parts by the dashed line in 
Fig. [7J the northeastern (NE) and southwestern (SW) regions. The NE region includes both 
the eastern and central clumps and the SW region includes the western clump and its south 
extension. The SW region only contains 2 or 3 IR galaxies, while the NE region contains 
at least 9 IR galaxies. The ratio of the IR galaxies to the number of cluster candidates 
selected by photo-z in the SW and NE region is 2 ± 1% (3 vs. 125) and 7 ± 2% (9 vs. 130) 
respectively, compared to 10 ± 2% (22 vs. 213) in the outskirt region (90" away from the 
cluster center). Although it may arise partly from projection effects, the difference between 
these two regions is statistically significant. In addition, four out of the five brightest IR 
galaxies of the whole cluster correlate to the NE region. These observations probably suggest 
a recent star-forming episode in the NE region. 

The X-ray study of the cluster may provide some support to this sc enario. Unlike the 
other two maj or clumps, the eastern clump is absent from the X-ray map (IClowe et al.ll2000l ; 



Jeltema et al 



200ll ; |Jee et al.ll2005al ). Based on this fact and th e anomalo u s X-ray profile of 



the central peak, as well as the temperature map of the region, |Jee et al.l (j2005al ) proposed 
that the eastern clump has passed through the dense region of the central clump recently 
as an off-center collision, during which the intracluster gas was stripped. It is possible that 
this recent off-center collision between the eastern and central clumps triggered SF in some 
of the IR galaxies. 



By comparison, the merging between the western and central clumps is pr obably at 
a very late stage, suggested by the lack of a shock-heated region between them (IJee et al. 



-20- 



2005al ). We define a central 9" x 25" rectangle as the interface of these two major clumps 
(solid box in Fig. [7]); we found two out of 21 photometrically selected cluster candidates 
with detectable IR emission in this region, a fraction of 10 ± 7%, comparable to the fraction 
in the outskirt region of the cluster. However, we found no spectroscopically confirmed 
IR member in this region. The uncertainties of the photometric cluster candidates and the 
projection effec t s furth er weaken the evidence for activity at this interface. As a comparison, 
Marcillac et al.l (120071 ) found no evidence of merging triggered SF activities in the interface 
of the two main merging clumps in the cluster RXJ0152 at a similar redshift. Combining 
the two results, clump interfaces do not appear in general to be sites of strongly enhanced 
star-forming activity. 



Tran et al.l (120031 ) studied E+A galaxies, characterized by strong Balmer absorption and 



little or no [O II] A3727 emission, in this cluster. The strong Balmer absorption is evidence of 
recent substantial star forming activities (< 2 Gyr) and the lack of [O II] emission indicates 
that those galaxies probab ly have no current SF, and therefore they are usually classified 
as post-starburst galaxies (IDressler fc Gunnlll983l ). There are 19 E+A cluster members in 
the region covered by IR observations and they are shown as the open stars in Fig. [71 Most 
of these galaxies have no IR emissi on and they are real po st-starburst galaxies in which SF 
has ceased at least a few Myr ago (jPoggianti & Wull2000l ). One third of the post-starburst 
galaxies are concentrated in the NE region of the cluster, three or four are related to the 
SW region, and the rest are scattered outside of the main body of the cluster. For the few 
post-starburst galaxies related with the SW region, only one is located in the very inner 
region, and the rest are distributed along the boundary region. This result suggests the NE 
region of the cluster has been active in SF for many hundreds of Myr, and the SW region, on 
the contrary, has been quiescent for a long time. It is consistent with the scenario suggested 
by X-ray analysis, that the central and eastern clumps experienced a recent collision, while 
the merging between western and central clumps is at a late stage. 

There are three E+A galaxies with IR emission, indicating ongoing SF. They are not 
post-starburst but rather dust-enshrouded starburst galaxies in which young stars are heavily 
obscured and their emission lines are extincted. The strong Balmer absorption suggests there 
are also a large number of A type stars, that have probably already moved out of the dusty 
star- forming regions, which indicates the SF in these galaxies started at least a few Myr 
ago (jPoggianti fc Wull2000l ). The dusty star-forming galaxies are all distributed outside of 
the main body of the cluster where gas stripping is not effective and the gas fuel for SF in 
the galaxies is not strongly depleted. This distribution helps explain why the SF in these 
galaxies continues for so long. 
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4.5. Mergers and Morphologically Irregular Members 



One remarkable feature of MS 1054-03 is its high fraction of merging galaxies. Van 
Dokkum et al. (1999, 2000) reported 13 ongoing mergers in this cluster, comprising 17% 
of the cluster population with L > L*. They also classified one more galaxy (H1532) as a 
merging candidate but did not include it in their merger list due to the lack of spectroscopic 
data for its companion. The photometric redshift suggests its companion is also a cluster 
member, and therefo re we inclu d e it a s a merger; this brings us to a total of 14 merging 
systems. In addition, iTran et al.l (120051 ) found 10 bound red galaxy pairs (5 of them already 
in van Dokkum's sample) with projected distance smaller than 30 h^ 1 kpc and relative 
velocity 8v < 300 km s" 1 . Most of the mergers between field galaxies a re accompanied by 



triggered SF (e.g., iLiu &: Kennicuttlll995l ; iPatton et al.l 120051 ) . However, ITran et al.l (120051 ) 



found that most of the merging galaxies in MS 1054-03 have no detectable [O II] emission 
lines and have probably lost their gas long ago. 

As pointed out in the previous section, [O II] emission may underestimate the SFR in 
galaxies due to dust extinction and is a less robust indicator of SFR than the IR luminos- 
ity. Therefore, we discuss the IR properties of these merging galaxies. Among 14 merging 
systems, 4 have 24 /im emission / 24 > 50 /iJy, H4683.4741, H6567, H2710 and H1532 (the 
last merger is the one confirmed by the photometric redshift). The brightest one is a double 
nucleus, highly disturbed disk system. The other three are in teracting pairs. One galaxy 
(H4822) common to two red pairs found by ITran et al.l (120051 ) has f'24 > 50 /xJy. There is 
another galaxy in those red pairs having weak 24 /im emission (^24 ~ 40 /xJy, not included 
in our spectroscopic sample). Altogether, about 29 ± 16% of the merging systems show 
detectable IR emission. For red pairs, this fraction is even lower, only about 10-20%. 



Postman et al.l (120051 ) classified the morphological types of spectroscopically confirmed 



cluster members of MS 1054-03. We correlate the IR galaxies in the spectroscopic sample 
with their classification. Fig. [8] shows the ACS images of these galaxies. Without distinguish- 
ing if galaxies are in merging systems, we found 21 ± 12% (4), 63±23% (12) and 16±10% 
(3) of the IR galaxies ar e early, late and i rregul ar type galaxies, respectively. These fractions 
are very similar to what iMarcillac et al.l (120071 ) found for the IR galaxies in RXJ0152. If we 
consider all the members of MS 1054-03 with morphological classification, 4±2%, 52±19% 
and 50±35% of early, late and irregular type cluster gal axies have detec t able IR emission 
(^24 > 50 /iJy). Again, these results agree with those of IMarcillac et al.l (120071 ) within the 
statistical errors. Therefore, only a very small fraction of early type galaxies in the spec- 
troscopic sample have strong SF, but half of late type and irregular galaxies have SFRs 



> 10 M t 



yr 



In addition to those galaxies in merging systems, galaxies with irregular morphologies 



-22 - 



may also be considered as experiencing interactions. Therefore, if we count the irregular 
galaxy (H4389) as a merger too, 6 IR galaxies out of 19 (32±15%) are related to galaxy 
merging/interaction. Furthermore, two IR galaxies (H6065 and H6372) classified as normal 
late type galaxies also show some irregular features. If their irregular features are also 
related to galaxy merging/interaction, then 42±18% of the IR galaxies may have galaxy 
merging/interaction triggered SFs. 



5. Discussion 



5.1. Comparison with RXJ0152. 7-1357 

From §4.2, we found that the IR LF of MS 1054-03 evolves as strongly as the field IR 
LF and has an over abundance of IR galaxies down to logLjR = 44.6, comparable to the field 
IR LF at the same redshift. This result is quite surprising given studies showing decreased 
SF in rich clusters. However, with only one cluster, it is not clear if the result is typical of 
rich clusters at z ~ 0.8 or peculiar to MS 1054-03. 



The study of IR galaxies in RXJ0152 by iMarcillac et al.l (120071 1 provides a good com- 
parison. The 24 /iin data for these two clusters have similar sensitivities and the optical data 
cover both central 5' x 5' regions of the clusters. X -ray data show that RXJ0152 also has two 
regions with peak emission (e.g.. iHuo et al . 2004) , and the dyna mical analysis of the cluster 
suggests an ongoing merger in the system (iGirardi et al.l 120051) . In additio n, photometric 



Tanaka et al. 



(2006) discovered two 



Marcillac et all (120071 ) 



and spectroscopic surveys by lKodama et al.l (120051 ) and 
large-scale filament-like structures hosting the central main cluster, 
found 22 IR galaxies in RXJ0152 confirmed by spectroscopic data as cluster members. We 
plot the IR LF deduced from their data as the open stars in Fig. [5j We did not make 
any incompleteness correction for the IR LF of RXJ0152 because their sample only includes 
galaxies with / 2 4 above the 80% completene ss limit and the spectr oscopic data in this cluster 
are quite deep, complete down to R = 24 (IDemarco et al.ll2005l ). Therefore, the IR LF of 
RXJ0152 without any correction is comparable to that of MS 1054-03 with incompleteness 
corrections above the 80% completeness limit at 24 /xm. 



As stated in §3.5, IMarcillac et al.l (120071 ) use a different set of SEDs to deduce a total 
IR luminosity from the 24 /im flux. Although it does not result in a significant difference, 
it is better if we compare the IR LF of RXJ0152 with that of MS 1054-03 deduced from the 
same method (open and filled triangles in Fig. [5]). There is a large difference in the faintest 
data points of the IR LF of RXJ0152 and MS 1054-03, due to the different cutoffs of 24 
jum flux density in the two samples: the RXJ0152 sample is only determined down to the 
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80% limit (^24 > 83 /iJy) while the MS 1054-03 sample, after incompleteness correction, is 
derived down to the 50% limit (/24 > 50 /iJy). The dashed vertical line in Fig. [5] shows 
the IR luminosity corresponding to the 80% limit. The lowest bin is partly below this 
limit. For both of the two brighter data points, RXJ0152 has a higher value than MS 1054. 
However, the differences are still within one sigma Poisson errors. The brightest data point 
for RXJ0152 includes 5 IR galaxies, while the data point for MS 1054-03 only includes 2 IR 
galaxies from the spectroscopic sample. The photometric selection adds two more galaxies 
in MS 1054-03 within this luminosity range, but as discussed before, those two are probably 
background sources. The IR LF of RXJ0152 confirms the strong evolution trend and the 
over abundance of bright IR galaxies we found in MS 1054-03 and it shows an even larger 
number of the brightest IR galaxies as pointed out by lMarcillac et al.l (120071 ). 



Marcillac et al.l (120071 ) find that most of the IR galaxies of RXJ0152 are distributed 



outside of the two major clumps indicated by X-ray emission. They also find a larger 
median redshift for these galaxies compared to the cluster redshift, which is identical to the 



l arger median redshift of the infalling late-type cluster members found by iBlakeslee et al. 
( 120061 ). Based on these facts, they suggest that infall of galaxies is probably responsible 



for much of the star formation activity we see in the cluster. In MS 1054-03, we did not 
find a difference in the redshifts of the IR galaxies compared with the rest of the cluster 
members. However, about 60% of the IR galaxies are located outside of the main body 
of the cluster (major clumps E, C and W). The projection effect may make the fraction 
even larger. Some of these galaxies probably correspond to the infalling galaxies found by 



Marcillac et al.l (120071 ). As discussed in §4.3 and §4.4, it is unlikely that the cluster only 
passively accretes star forming galaxies from the surrounding field and those galaxies have 
a high level of SF simply due to their recent origination from the field. On the contrary, 
it is very possible that we are seeing an increased SFR in infalling galaxies triggered by 
the galaxy-IGM interaction as shown by prev i ous t heoretical and observational evidence 



(IFuiita & Nagashimalll999l : iGavazzi et al.lll995l . |200l|). 



However, we can not rule out other mechanisms being responsible for the SF in these 
galaxies. In MS 1054-03, about one third of the outside IR galaxies are associated with 
some minor clumps of the cluster. Their SF m ay arise from pro cesses more common in 
the group environment, e.g., galaxy interactions (ILewis et al.ll2002l ). For the 40% of the IR 
galaxies correlated with the main body of the cluster, the majority are associated with the 
NE region of the cluster where a collision of subclumps might have occurred recently. This 
result indicates the interaction of the subclumps and the processes accompanied with it, e.g., 
the tidal gravit ational field , may also play a role in triggering SF and cause the concentration 
of IR galaxies (iBekkil Il999h . 
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The difference in the mass of the two clusters can com plicate the compa rison. Using 
the same weak-lensing technique and data of similar quality, IJee et al.l (j2005ai lbl) produced 
enclosed mass profiles for both of them. The profiles show that MS 1054-03 is much richer 
and more massive than RXJ0152. The enclosed mass within 1 Mpc of MS 1054-03 is about 
two times as large as that in RXJ0152. Therefore, with the similar 5' x 5' IR and spectro- 
scopic/photometric coverage of the clusters, we actually only observe the inner part of MS 
1054-03 but reach the outside region in RXJ0152, where most of its LIRGs reside. It is pos- 
sible that we would find more infalling IR galaxies if our IR and spectroscopic/photometric 
data extended further to the outer regions of MS 1054-03. The many IR galaxies distributed 
at the very edge of the survey region (see Fig. [7]) seem to support this argument. This would 
also explain the slightly higher IR LF of RXJ0152 compared with MS 1054-03. 

Although both MS 1054-03 and RXJ0152 have two X-ray peaks and two corresponding 
major clumps indicating a me rger, MS 1054-03 lacks shock-heated regions between the two 
X-ray peaks (IJee et al.l l2005al ). bu t RXJ0152 has exces s X-ray emission between the two 
clumps suggestive of a shock front (IMaughan et al.ll2003l ). The differences may indicate the 
different merger stage the two clusters are in: MS 1054-03 is probably at a post-merger stage 
while RXJ0152 is at a pre-merger stage. This difference may contribute to the slightly larger 
star formation rate in RXJ0152 than in MS 1054-03. 



5.2. Evolution of the Integrated SFR in Clusters 



We have already compared the IR LF of MS 1054-03 to the Coma cluster and found 
a strong evolution in both 0* and L* IR . Another way to compare the SFR in clusters is 
to compare the ir integrate d SFR s with in a certain radi u s. Us ing th e SFR measured from 
emission lines, iFinn et al.l (120041 . 120051 ). iKodama et al.l (120041 ). and iHomeier et al.l (120051 ) 
compared the integrated SFRs of several clusters within 0.5i?2oo- They also compared the 
integrated SFRs normalized by the cluster masses. The mass-normalized integrated SFRs 
are comparable to the fractions of star-forming galaxies in the clusters, which are widely 
used in many systematic studies of cluster SF. Some suggestive correlations between the 
integrated SFRs and redshifts/masses of the clusters were found. However, the results are 
very uncertain. 



There are concerns about the usu al methods for e s timat ing the cluster masses. IFinn et al 



( 120051 ) used velocity dispersion, while IHomeier et al.l (120051 ) suggest X-ray temperature may 
be a better indicator of the mass. However, velocity dispersion and X-ray temperature are 
valid mass estimators only in relaxed clusters under the assumption of hydrodynamic equi- 
librium, which is often questionable, especially for clusters at high redshift. To clarify the 
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results found in those studies, we add more data points by including the integrated SFRs of 
three clusters observed by MIPS (Coma, MS 1054-03, RXJ0152), and four clusters observed 
by ISOCAM (A2218, A1689, A2219, CI 0024). The seven clusters with the SFRs measured 
from Ha emission are also added (A1367, AC114, A2390, CI 0023, CI 1040, CI 1054 and CI 
1216). Despite the systematic difference between the SFR measured from the emission line 
strengths and from the IR luminosity, there is general agreement after extinction correction. 
We also add an average value for clusters at 0.3 < z < 0.5 deduced from IRAS data. For 
the cluster mass, we take the mass measured from lensing analysis whenever it is available 
because it is free from any assumptions about the dynamical state of the clusters. We limit 
our calculation to within the 0.5i?2oo region and set the cutoff in the SFR as 2 M yr _1 . 
The references and details about the integrated SFR and the masses of these clusters are 
provided in Appendix A. 

In panels a and b of Figj9j we plot the integrated SFRs as a function of redshift and clus- 
ter mass. In panels c and d, we plot the mass-normalized values. The integrated SFRs show 
a weak evolution with redshift. The evolution is more pronounced in the mass- normalized 
SFRs, approximately oc (1 + z) 5 . The nonparametric Spearman tests show that the signifi- 
cances of the correlations are 97% and 99% for the integrated SFRs and the mass-normalized 
integrated SFRs. However, this evolution trend is complicated by the anticorrelation between 
the mass-normalized integrated SFRs and the cluster masses. Although the integrated SFR 
does not show an apparent correlation with mass, the mass-normalized one has an anticorre- 
lation with mass, oc M~ ' 9 , with a Spearman sign ificance of 97%. This anticorrelation ag rees 



with the results found in previous compariso ns (IHomeier et al. 12005c iFinn et al.l 120051 ). It 



also agrees with the anticorrelation found by Poggianti et al. ( 20061 ) between the fraction 



of star-forming galaxies and the cluster velocity dispersion, in the sense that the fraction 
of star-forming galaxies is comparable to the mass-normalized integrated SFR. Given this 
anticorrelation, the evolution we found in the mass-normalized integrated SFR is probably 
due to the different masses of the low-redshift and high-redshift clusters in our sample. 

As shown in panel e of Figj9l the clusters with z < 0.5 in our sample are on average more 
massive than those with z > 0.5. Therefore, the increased mass-normalized integrated SFRs 
at higher redshifts found in the sample could merely be a selection effect. To disentangle the 
mass factor from the evolution trend, we select a subsample of clusters with a mass range of 
3 x 10 14 M Q <M<12x10 14 M q , in which both low-redshift and high-redshift clusters have 
good sampling, and plot their mass- normalized SFRs vs. redshift in panel /. The evolution 
trend of the subsample becomes much weaker, but it still has a significance of about 89%. In 
addition, even if the evolution of the mass-normalized integrated SFR can be largely caused 
by the anticorrelation between the mass-normalized integrated SFR and the cluster mass, 
the evolution of the integrated SFR without mass normalization is still significant and can 
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not be easily explained by the mass selection effect. 

Among all the clusters in our sample, CI 0024 has the largest integrated SFR with an 
intermediate redshift and cluster mass. Its integrated SFR is at least five times larger than 



those clusters with similar masses. Using MIPS 24 /zm data, iGeach et al.l (120061 ) also found 
a very significant excess of mid-infrared sources up to r < 5 Mpc in CI 0024 compared to 
another cluster MS 0451-03 at z = 0.55. These results suggest that CI 0024 is quite unusual 
compared to other clusters in the sample. Another unusual aspect of CI 0024 is its relatively 
faint X-ray emission compared to its large mass. It has a mass a little larger than RXJ0152, 
but its X-ray luminosity is only about a fifth of RXJ0152. 

The eastern clump of MS 1054-03 also has significant star formatio n activities but i s 



absent in X-ray emission, contradicting the expected Lx from its mass (IJee et al.ll2005al ). 
For a comparison, we calculate the integrated SFRs for the eastern clump of MS 1054-03 
separately and plot it as an open star in Fig.[9j This result has the interesting implication that 
clusters or clus ter subclumps with u nusually low X-ray emission may h ave very active SF. A 



recent work by lPopesso et al.l (120071 ) seems to support this conclusion. iPopesso et al.l (120071 ) 
studied 137 Abell clusters and found that clusters with lower X-ray luminosity than expected 
from the Lx — M relation, the so-called X-ray-Underluminous Abell clusters (AXU), show a 
velocity distribution characteristic of accretion and have a higher fraction of blue galaxies in 
their outer regions. They suggest the low X-ray luminosities of these clusters are due to the 
ongoing accretion or merging process. Although CI 0024 is not exactly X-ray underluminous 
according to their definition, the exceptionally high SFRs in it and in the eastern clump of 
MS 1054-03 generally agree with the scenario they propose. They also found that about 
40% of the clusters they studied are AXU, indicating that AXU clusters are not a small 
minority, at least at z < 0.4, and suggesting they probably host more SF in total than 
the X-ray-luminous clusters. Since most of the clusters in our sample are X-ray luminous 
ones, the SFRs have probably been biased towards the lower value and have an evolution 
reflecting only conditions in well relaxed systems with substantial amounts of hot, X-ray 
emitting plasma. 

As pointed out previously, the star-forming galaxies tend to be located in the outer 
regions beyond 0.5_R 2 oo- Thus, the integrated SFRs within 0.5-R 2 oo only present a portion of 
the total star forming activities in the cluster. This cutoff effect is especially significant for 
RXJ0152 due to its irregular morphology. The integrated SFR within 0.5i?2oo of RXJ0152 
only accounts for about 13±3% of its total SFR in the survey region, much smaller than 
the 50% expected from a singular isothermal (SIS) distribution of the star-forming galaxies. 
For MS 1054-03, this fraction is higher, 70%, comparable with 66% expected from a SIS 
distribution. 
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5.3. Ram Pressure Stripping 



In the preceding section, we showed that the mass-normalized integrated SFRs of clus- 
ters have an anticorrelation with the masses. This result indicates that massive clusters are 



prob ably more effective in suppressing SF than the low mass ones (see also, iPoggianti et al. 



20061 ). We also found that very few IR galaxies in MS 1054-03 are distributed in the region 
of high mass density, especi ally in the southwes tern part. On the other hand, gas stripping 
by ram pressure in clusters (IGunn fc Gottlll972l ) is also found to be mo re pronounced in the 
mass ive clusters and more effective in the high density regions (e.g., iGiovanelli fc Haynes 
19851 ). The coincidence may suggest SFR suppression in the clusters due to gas stripping by 
the ram pressure of the intracluster medium (ICM). 



T o investigate the effect of the ram pressure gas stripping, followin g iHomeier et al. 
( 2005 ). we calculate the effective radius for this process. Jee et al. (|2005a ) found the X-ray 
surface brightness profile of the cluster is best fitted by an isothermal (3 gas model with 
(3 = 0.78 ± 0.08 and r c = 16" ± 15". Due to the complex morphology of the cluster, this 
isothermal form does not fit the surface brightness profile of the inner r < 45" region, but 
it fits the outer region very well and predicts a projected mass profile consistent with the 
result from the weak-lensing analysis. The virial ra dius of the cluste r is 1.7 ± 0.2 Mpc and 
the corresponding virial mass is 1.2 ± 0.2 x 10 15 M (I Jee et al.ll2005al). If we assum e the gas 
mass is about 10% of the total mass, as suggested in iNeumann fc Arnaudl (120001 ) . we can 
deduce a central gas density po, gas = 8 x 10 14 Mpc~ 3 . Acco rding to the ram pressure 
stripping criterion for the gas in a disk (IFujita &: Nagashimalll9991 ) and the gas profile of MS 
1054-03, we obtain a ram pressure stripping effective radius of r rp = 63" (~ 0.5 Mpc) for a 
Milky Way-type galaxy with a velocity of 1000 km s -1 . 



We show the r rp in Fig. For RXJ0152, r rp ~ 0.3 Mpc (IHomeier et al.ll2005l ) and none 
of its star- forming galaxies lies within this radius (also see the discussion in Marcillac et al., 
2007). In MS 1054-03, the majority (80%) of the IR cluster galaxies are distributed outside 



of r rp , providing strong evidence of SF suppression due to ram pressure stripping. For the 
six IR galaxies within r rp , only two of them are brighter than / 2 4 — 80 /iJy. The fraction 
of IR galaxies (4 ± 2%) in this region is about half of that in the outskirt region but still 
larger than the fraction in the southwestern part, which supports the suppression effect of 
ram pressure stripping, but also suggests it is less effective for some galaxies and r rp is only 
an approximate measurement of its effectiveness. There may be a number of explanations. 
First, this could just be a projection effect, that is at least some of the six galaxies lie in 
front of or behind the ram stripping region. Second, the isothermal density profile is only 
an approximate description of the gas distribution, so r rp is only a rough indicator of the 
effectiveness of ram pressure stripping. Third, the ram pressure stripping criterion we used 
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to deduce r rp is only for the gas in the disk of a galaxy. The gas in the inner disk of a 
galaxy is harder to strip and therefore any SF occurring there, e.g., circumnuclear SF, is 
more difficult to suppress. 

Another interesting fact is that about half of the E+A galaxies are distributed along 
the effective radius. This behavior suggests gas stripping as the reason that the star forming 
activity stopped in these galaxies. 



5.4. Galaxy Interaction 

Many studies find that most of the star- forming galaxi es in both the field and clusters 



are m orphologically not strongly disturbed at z ~ 0.8 (e.g., iBell et al.ll2005l ; iMarcillac et al. 



20071 ) and argue that mechanisms that would dramatically disturb the morphology of a 
galaxy, e.g. strong galaxy interactions, can not be the major factor causing the change of 
SFRs with epoch and environment. However, in MS 1054-03, we found 5 out of 19 IR galaxies 
of our spectroscopic sample are in interacting systems, one isolated galaxy is irregular, and 
two more have some irregular features. Therefore, more than 30% of the star-forming galaxies 
in the spectroscopic sample may be related to galaxy interactions. Among 5 IR galaxies in 
interacting systems, only two have disturbed morphology, and the other three look regular. 
If we only count the irregular galaxies as mergers, no matter if they are in merging systems 
or not, only 16% (3 out 19) of the IR galaxies are related to galaxy interaction, consistent 



with the result from IBell et al.l (120051 ). who found less than 30% of field IR galaxies are 



strongly interacting at z ~ 0.7. 

For the regular IR galaxies in interacting systems, it seems that the galaxy interac- 
tions may have triggered their star formation activity but do not change their morphology 
sufficiently to be obvious. This could be due to the different time scales on which SF and 
morphological change occur during a galaxy interaction. It is also possible that those interac- 
tions are only strong enough to trigger the instability of the galaxies and cause star formation 
activity, but not to cause observable morphology distortions. Interaction triggered SF is not 
always accompanied by disturbed morphology and the study of galaxy morphology at the 
redshift of MS 1054-03 is probably only sensitive to the strongly interacting systems. 

On the other hand, the majority of the bright interacting systems (~ 70%) in the 
cluster do not have detectable IR emission. Some of them have strong interacting features, 
e.g., double nuclei, distorted morphologies, but these characteristics are not accompanied 
by strong SF. This is probably because most of the interacting galaxies have already lost 
their gas (dry merger) while falling into the cluster due to, e.g., ram pressure stripping, and 



cannot support a high level of SF. An example of a d ry merger with litt le star formation 
increase has also been found among local galaxies (e.g., iBoselli et al.l 120051 ). 



6. Conclusion 

Using the MIPS 24 /zm data for the rich cluster MS 1054-03 at z = 0.83, we found 
19 IR emitting cluster members selected by spectroscopic data and 15 additional IR cluster 
member candidates selected by photometric data. 

We constructed the IR luminosity function of the cluster and find a strong evolution 
when compared with the IR LF of the Coma cluster with a similar mass at z = 0.02. The 
characteristic IR luminosity {L* IR ) of MS 1054-03 is about one order of magnitude larger than 
that of the Coma cluster. The SFR density integrated from the IR LF is about 16 times 
larger than that in the Coma cluster. The evolutionary trend of the IR LFs from Coma to 
MS 1054-03 is similar to the evolution of the IR LFs in the field. The comparison of the 
mass normalized integrated SFR of MS 1054-03 with several other clusters seems to agree 
with the evolution suggested by the IR LFs, but it is less conclusive because of the combined 
mass and redshift dependence of the SFR. The similar SFR evolution in the clusters and in 
the field favors some internal mechanism, e.g., the consumption of the gas fuel in galaxies, 
as being responsible for the decline of SFR in different environments. 

A substantial fraction (13±3%) of cluster galaxies are forming stars actively. Although 
the fraction is lower than that in the field (52%), the overdensity of the IR galaxies in the 
cluster is still quite high, ~20. Such a high level of SF is evidence against the scenario 
that the cluster is only accreting star-forming galaxies from the surrounding field passively, 
after which their star formation is quenched. Instead, it appears that many cluster galaxies 
continue to form stars at a high rate. A number of cluster galaxies still have large amounts 
of gas and their SF can be triggered by the interactions with the intergalactic medium, 
with other galaxies, or, by tides. However, there are few IR galaxies distributed in the high 
density regions of the cluster, indicating the suppression effect of ram pressure stripping on 
the SFR in those regions. Both the IR galaxies and the E+A galaxies of the cluster show 
a concentration in the NE region of the cluster, supporting the scenario that an interaction 
between subclumps occurred recently and enhanced the SFR. 

About half of the bright late type and irregular cluster galaxies have detectable IR 
emission, but for early type galaxies this fraction is only about 4%. Only 29% of the mergers 
in the cluster have detectable IR emission. The majority of the mergers probably have lost 
their gas fuel long ago and can not support a high level of SF. More than 30% of IR galaxies 
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show evidence of galaxy interaction, and only half of them have irregular morphologies, 
suggesting the interaction-triggered morphological change and star formation activities of 
galaxies have different time scales and intensities. 
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A. Integrated SFRs and Masses 



The integrated SFR of MS 1054-03, RXJ0152 and Coma, are 372, 134 and 35 M yr" 1 , 
respectively. Since the MIPS IR data are only complete down to ~ 10 M yr" 1 for MS 1054- 
03 and RXJ0152, we have applied a correction factor of 1.5 estimated from their best fitted 
IR LFs to the obs erved integrated SFRs. We adopt weak-lens ing masses of 1.1±0.1 x 10 15 M Q 
Jjee et al.ll2005ah and 4.5 ± 2.7 x 10 14 M Jjee et alJl2005bh for MS 1054-03 and RXJ0152. 
For the Coma cluster, we use a mass of 1.4 ± 0.4 x 10 15 M Q from the dynamical analysis 
Jlokas fc Mamonll2003h . 



The four c lusters observed by ISOCAM are A2218 flBiviano et al.ll2004f) . A1689 flDuc et al 



2002|), A2219 flCoia et al ll2005bh an d CI 0024 flCoia et al.l l2005al) We adopt their masses 
as 4.8 ± 1.4 x 1Q 14 JPratt et al.l liooih. 1.93 ± 0.2 x 1 15 faroadhurst et aD l2005h . 
1.0 ± 0.7 x 10 15 M JPahle!l2006h and 5.7 ± 1.1 x 10 14 M fcneib et al.ll2003h . These are all 
lensing masses, except for A2218, which is measured by fitting the X-ray temperature pro- 
file. The ISOCAM surveys of A2218 A1689 and A2219 only covered the central ~ O.2.R200 
regions. Following iFinn et al.l (120041 ). we correct for the small coverage by assuming the 
star-forming galaxy distribution has a singular isothermal (SIS) profile. This correction will 
give a lower limit to the integrated SFR if the star-forming galaxies are more likely to reside 
in the outskirt region of the cluster as suggested by MS 1054-03 and RXJ0152. For A2219 
and CI 0024, we also need to correct for the incompleteness of the detection limits, which are 
both ~ 10 M & yr _1 . Since we do not have IR LFs for these two clusters, we use the IR LFs 
of field galaxies at similar redshifts to estimate the correction to extend the detection limit 
down to 2 M Q yr^ 1 . This correction is made on the assumption that the shape of the IR LF 
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of the rich cluster does not differ significantly from that of field galaxies at the same redshift, 
which is probably true given Coma and MS 1054-03 as examples. After these corrections, 
the integrated SFRs for A2218, A1689, A2219 and CI 0024 are: 14, 64 ± 17, 307 and 753 M 
yr -1 . For the SFRs deduced from the IR luminosity, we assume a 50% error if it has not 
been given, considering the typical uncertainties of the IR flux measurement. 

We include seven clusters with SFRs me asured from Ha em ission in the comparisons, 
A136 7 Jlglesias-Paramo et al.ll2002h . AC114 JCouch et alJbpOlh. A23 9 faalogh fc Morris 
2000h . CI 0023 JFinn et all booT and three clusters from iFinn et~aD J2OQ5J): CI 1040, CI 
1054, CI 1216. The SFRs of these clusters are measured from Ha narrowband imaging, 
except AC114 which is measured from Ha spectro scopy. The mass o f A1367 is estimated 
using the dynamical analys is, 7.1 ± 1.5 x 10* 4 M ( IGirardi et al.lll998l ). We use the lensing 
mass of 7.3+ 4 ; 4 x 10 14 M rtNataraian et al.lll998l) for AC1 14. For A2390, we adopt a mass 
of 13.6 ± 0.7 x 10 14 M Q from X-ray analysis (lAUen et al.ll200ll ). whic h gives a consisten t 
result wit h the lensing analysis in the inner region where it is available (ISquires et al.lll996l ). 
Follo wing iKodama et al.l (120041) . we use a dynamical mass of 2.3 ±1.2 x 10 14 M & for CI 
0023 jPostman et alJll99«h . For CI 1040, CI 1054, and CI 12 16, we use the len sing masses of 
0.55i°vg x 10 14 M , 4.8±i;| x 10 14 M and 9.5+}| x 10 14 M Jdowe et alil2006h . respectively. 
For A13 67, we correct for the incom plete coverage both in space and velocity range by a factor 
of 1.7 (jlglesias-Paramo et al.l 120021 ) and obtain an integrated SFR of 29.1 ± 3.1 M yr -1 . 
For AC114, we also include a correction fa ctor of 2.8 for t h e sam pling fraction and aperture 
bias of the spectroscopic survey following IKodama et al.l (120041 ) . All of these Ha surveys 
are complete down to 2 M yr -1 and we make a correction for incomplete survey coverage 
for CI 1216. The integrated SFRs for AC114, A2390 and CI 0023 are 21.6 ± 19.5, 80 ± 28 
and 71 ± 23 M yr -1 . These values are slightly smaller than those given by iHomeier et al. 



(120051 ) because we limit the integration to galaxies with SFR > 2 M yr 1 . For CI 1040, CI 
1054 and CI 1216, we obtain integrated SFRs of 63 ± 12, 90 ±19 and 369 ±55 M. yr" 1 . The 
results for CI 1040 and CI 1216 are similar to those found by IFinn et al.l (120051 ) . However, 
because the lensing mass o 
the velocity dispersion by 



larger than that given by IFinn et al 



C11054 is a bout three times larger than the mass estimated from 
(2005|), the integrated SFR is also about three times 
(120051 ) due to the larger i^oo- 



Finn et a 



In addition, the SFRs of these seven clusters are deduced by assuming a Ha extinction 
of 1 mag, which corresponds to Ay ~ 1.2 mag. However, as mentioned in section 4.3, most 
of the star-forming galaxies have extinctions larger than this value. For example, all the IR 
galaxies in MS 1054-03 have Ay 1 larger than 1.5 mag. To account for the underestimate 
of the extinction for the Ha deduced SFRs, we need to know the average Ay of galaxies in 
each cluster. Since we do not have the Ay measured for these galaxies, we have to assume 
the shapes of the IR LFs of clusters are not very different from those of field galaxies at 
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the same redshifts and estimate the average IR luminosities for those galaxies with SFR 
> 2 M Q yr -1 . This assumption may be inaccurate, but both Coma and MS 1054-03 IR LFs 
seem to support it. From the average IR luminosit ies, we can estima te the A\P with the IR- 



luminosity-dependent extinction formula given by lChoi et al.l (120061 ). Comparing those Ay 1 
with the assumed 1 mag extinction at Ha, we deduce correction factors for the integrated 
SFRs in these clusters, which range from 1.2 to 1.5. We apply the corrections for all the Ha 
deduced integrated SFRs when comparing them with those deduced from IR luminosities. 



Using IRAS 60 /mi data, iKelly &: Riekd (119901 ) find the average flux of a sample of 



clusters by stacking and folding their addscan signals. There are 58 clusters in their sample 
with 0.3 < z < 0.5, and their average flux is about 29.1 mJy (They also calculated a 
value for a sample of local clusters with z ~ 0.05, but due to the uncertainties caused by 
the different techniques used to deduce the average flux, we do not consider those clusters 
here.) Each scan of the z > 0.3 clusters is about 20' long and 5' wide and the signal is 
recorded as the one dimensional flux along the scan. The width of the scan is approximately 
equal to -R200 of a cluster with 7 x 10 14 M & at this redshift and the PSF fitting of the scan 
signal along the length includes the flux approximately from the central 5' region. So the 
average flux they obtained comes from a similar region as our integrated SFR for a cluster 
with 7 x 10 14 M . For the clusters with larger/smaller mass than 7 x 10 14 M & , the flux is 
an underestimate/overestimate of the total flux within 0.5i?2oo regions due to the limitation 
of the scan region. Because the evolved Coma IR LF suggests that about one third of the 
total IR luminosity is coming from LIRGs at this redshift, we deduce a composite Lm/L eo 
ratio for the clusters by averaging the ratio given by LIRG and late type SEDs weighted 
by this factor. Using this composite Ljr/Lqo ratio, we can correlate the 60 /im flux of 
the cluster to the total IR luminosity and therefore the total SFR. The error of this total 
SFR is dominated by the uncertainties caused by the composite L IR /L m ratio. Since the 
average 60 fim flux of these clusters is obtained by stacking, it is not limited by the detection 
limit of the observations and the total SFR deduced from it also includes the contribution 
from the galaxies with SFR < 2M Q yr -1 . Those galaxies contribute 30% of the total SFR, 
estimated from the IR LF. Using this factor, we convert the average total SFR to the average 
integrated SFR. Assuming a typical mass of 7 x 10 14 M for these clusters, we obtain the 
mass-normalized integrated SFR. The assumption of the typical mass here will cause some 
uncertainty on this data point, but it will partly offset the effect caused by the limitation of 
the scan region. 
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Table 1. Cluster galaxies with 24 /im emission 



ID a 




f b 


logL IR 




(Mag) 


(/^Jy) 


(ergs s _1 ) 


Spectroscopic Sample 


1140 


20.24± 0.06 


363± 33 


45.37 


553 


19.83± 0.02 


317± 30 


45.13 


1663 


20.51± 0.05 


140± 28 


44.69 


791 


21.90± 0.09 


76± 44 


44.68 


102 


19.86± 0.06 


184± 32 


44.68 


909 


20.76± 0.01 


116± 24 


44.63 


1316 


20.91± 0.01 


96± 26 


44.52 


1200 


20.07± 0.07 


107± 36 


44.45 


581 


21.61± 0.03 


62± 30 


44.43 


1357 


20.86± 0.04 


82± 25 


44.32 


695 


20.85± 0.02 


76± 22 


44.29 


874 


20.68± 0.04 


68± 28 


44.25 


211 


19.97± 0.02 


63± 31 


44.22 


725 


21.50± 0.08 


51± 34 


44.21 


170 


21.14± 0.06 


59± 25 


44.19 


107 


20.04± 0.01 


56± 31 


44.17 


1108 


20.64± 0.04 


54± 25 


44.16 


166 


21.25± 0.01 


55± 24 


44.15 


195 


21.01± 0.05 


52± 17 


44.14 


Photometric Sample 


135 


21.19± 0.05 


510± 36 


45.59 


1098 


20.55± 0.04 


365± 33 


45.43 


870 


21.89± 0.06 


144± 15 


45.04 


1546 


22.09± 0.11 


89± 33 


44.82 


1528 


21.87± 0.17 


84± 39 


44.75 


1849 


20.85± 0.02 


124± 29 


44.72 


1121 


22.82± 0.26 


61± 21 


44.67 


788 


23.09± 0.21 


54± 10 


44.61 


828 


20.63± 0.04 


116± 25 


44.56 
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Table 1 — Continued 



ID a 




f b 


\ogL IR 




(Mag) 


OJy) 


(ergs s^ 1 ) 


950 


21.89± 0.10 


63± 23 


44.55 


420 


21.50± 0.13 


63± 27 


44.41 


278 


21.93± 0.11 


51± 34 


44.40 


1716 


20.38± 0.03 


88± 17 


44.37 


1750 


21.56± 0.11 


58± 29 


44.35 


1846 


20.62± 0.02 


76± 60 c 


44.30 



^Galaxy IDs and K^ar are taken from 



Forster Schreiber et al.l (120061 ) . 



b The 24 /im flux errors are estimated 
within a fixed aperture (r < 5") and are 
scaled by the ratio of the DAOPHOT PSF 
fitting flux vs. the fixed aperture flux. The 
method would overestimate the error for 
some faint sources. 

c The source has a low nominal SNR of 
1.3 due to the reason mentioned in note b; 
it has been visually confirmed as a secure 
detection. 
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Fig. 1. — The completeness of the IR and spectroscopic surveys. The left panel is the 
completeness at 24 /im. The sample is about 80% complete down to / 24 ~ 80 /zJy. The right 
panel is the completeness of the spectroscopic survey as a function of Igu magnitude. The 
dashed vertical line is approximately the 90% completeness limit of the photometric survey. 



-42- 




■1.5 -1.0 -0.5 

'og f Ks / f 24 



1.0 



Fig. 2. — The rest frame Lm/ L^/^i+z) ratio as a function of K s - 24 c olor. The open 



symbo ls are the data points deduced from the template SEDs given by iDevriendt et al. 



( 11999I ). Upward triangles, square, and stars denote the normal spirals, LIRG and ULIRGs. 
The filled circle and the downward triangle are the results of the interpolation of the galaxies 
of the spectroscopic and combined samples from their K s - 24 colors. The small dots are the 
Ljpi/ L 2 4/(i+ z ) ratio of galaxies deduced from the second method in §3.5. 
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Fig. 3. — The comparison of SFR/^ and SFRjo//]. The filled and open circles are the SFRs of 
the IR galaxies in MS 1054-03, with and without extinction correction for SFR[ //]. For the 
three IR galaxies without [O II] emission lines, SFR/^j are used as the upper limit of SFR[o//j 
and they are plotted as open circles with downward arrows. The filled and open triangles are 
the SFRs of the IR galaxies in RXJ0152, with and without extinction correction for SFRro/n. 
The dotted vertical line is the SFR approximately corresponding to the 80% completeness 
limit of the 24 fim observation. The solid line indicates the one-to-one correlation. 



-44- 




Fig. 4. — The comparison of SFR/^ and SFR^y. The open circles are the spectroscopi- 
cally confirmed IR cluster members, and the open squares are the ones selected by their 
photometric redshifts. The filled circles and squares are the results after applying extinc- 
tion correction. The dotted vertical line is the SFR approximately corresponding to the 
80% completeness limit of the 24 observation. The solid line indicates the one-to-one 
correlation. 



-45- 



!og L IR /L e 

10.0 10.5 11.0 11.5 12.0 

— i — i — i — i — i — | — i — i — -i — i — | — i — i — i — i — | — i — i — i — i — | — 

MS1054: q ^ Spec-z 




43.5 44.0 44.5 45.0 45.5 46.0 

log L, R [ergs s _1 ] 



Fig. 5. — The IR luminosity function of MS 1054-03. The open and filled circles are the 
result of the spectroscopic sample without and with spectroscopic and IR incompleteness 
correction. The open and filled triangles are the result of the spectroscopic sample deduced 
from the second method in §3.5, without and with incompleteness correction. The open and 
filled squares (shifted to the bright end by 0.04 for clarity) are the result of the combined 
sample without and with IR incompleteness correction. The open stars (shifted to the faint 
end by 0.04 for clarity) are the IR LF of RXJ0152 from iMarcillac et al.l fl2007l ). Since they 
only include galaxies with / 24 > 80 /iJy and the data are very incomplete at \ogLm < 44, 
we draw the faintest point as a lower limit. The solid curve is the best fitting Schechter 
function to the corrected spectroscopic IR LF. The dotted curve is the best fitting Schechter 
function to the IR LF of the Coma cluster. The dashed curve is the Coma LF evolved to 
z = 0.83 with the same evolution trend as the field IR LF. The vertical dashed line is the 
IR luminosity corresponding to the 80% detection limit at 24 //m. 
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Fig. 6. — The IR luminosity function of MS 1054-03 compared with the field IR luminosity 
function at z ~ 0.8. The filled symbols are from the spectroscopic data after incompleteness 
correction. The solid curve is the best fitting Schechter function. The dotted curve is the 
field IR LF at the cluster redshift. The shaded area presents its one sigma uncertainties. 
The vertical dashed line is the IR luminosity corresponding to the 80% detection limit at 24 
/zm. 




Fig. 7. — The spatial distribution of the IR galaxies in MS 1054-03. Open squares are the IR 
galaxies from the spectroscopic sample, and open triangles are the additional IR members 
selected by the photometric redshifts. The sizes of the symbols are proportional to their IR 



luminosities. LIRGs are also indicated by filled circles. Small open circ 
members selected by the photometric redshifts. E+A Galaxies select ed by 



es are the cluster 



Trail et all fl2003h 



are indicated as open stars. The contours are the mass contours from lJee et al.l (j2005al ). The 
eastern, central, western and four minor mass clumps are labeled as E, C, W, and M1-M4. 
The solid circle is the region with effective gas-stripping. The rectangle is defined as the 
interface region of the two major clumps. 
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Fig. 8. — The ACS images of spectroscopically confirmed IR galaxies in MS 1054-03. The 
size of each image is about 6" x 6". 
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Fig. 9. — (a) and (6), the integrated SFRs vs. redshifts and cluster masses; (c) and (d), 
the mass-normalized SFRs vs. redshifts and cluster masses; (e), the cluster masses vs. 
redshifts; (/), the mass-normalized SFRs vs. redshifts for clusters with 3 x 10 14 M < 
M < 12 x 10 14 Mq. The mass selection limits are indicated as the two dashed horizontal 
lines in (e). The filled stars are three clusters observed with MIPS: Coma, MS 1054-03, 
and RXJ0152. The open star is the eastern clump of MS 1054-03. Filled circles are the 
clusters observed with ISOCAM, and open circles are from Ha emission line measurements. 
The dotted curves in (c) and (/) are the fitted correlation between mass-normalized SFRs 
and redshifts for all the clusters (not including the eastern clump of MS 1054-03). The 
dotted curve in (d) is the fitted correlation between mass-normalized SFRs and masses (not 
including the eastern clump of MS 1054-03). 



